Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



CERN-PH-TH/2006-226 



CONSTRAINTS ON PARAMETERS OF RADIATIVELY DECAYING DARK MATTER FROM THE GALAXY CLUSTER 

IE 0657-56 

A. Boyarsky'-^'\ O. Ruchayskiy", and M. Markevitch^-* 

Draft version February 5, 2008 



o 
o 

(N 
> 

o 



ABSTRACT 

We derived constraints on parameters of a radiatively decaying warm dark matter particle, e.g., the mass and 
mixing angle for a sterile neutrino, using Chandra X-ray spectra of a galaxy cluster IE 0657-56 (the "bullet" 
cluster). The constraints are based on nondetection of the sterile neutrino decay emission line. This cluster 
exhibits spatial separation between the hot intergalactic gas and the dark matter, helping to disentangle their 
X-ray signals. It also has a very long X-ray observation and a total mass measured via gravitational lensing. 
This makes the resulting constraints on sterile neutrino complementary to earlier results that used different 
cluster mass estimates. Our limits are comparable to the best existing constraints. 

Subject headings: Dark matter — elementary particles — galaxies: clusters: individual (1E0657-56) — line: 
formation — neutrinos — X-rays: galaxies: clusters — X-rays: individual (IE 0657-56) 
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1 . STERILE NEUTRINO AS WARM DM CANDIDATES 

A number of works appeared recently on the sub- 
ject of a sterile (or right-handed) neutrino as a pos- 
sible dark matter (DM) c a ndida te (e.g.. lAsaka et all 

2006b'a; 
2006a; 




azauan 



lAsa ka et al. 
Kiemer-S0rensen et at] 



2005; Asaka & Shaposhniko 

Bovarskv et al. 2006a b.c 

Shaposhnikov & Tkachev) 1200^ 

2006; .Watson et al., ,20061). The discovery o f neutrino os- 
cillations (see, e.g., IStrumia & Vissanil |2006| for a review) 
made the existence of a sterile neutrino quite plausible and 
spurred the interest in this candidate. Several factors make 
dark matter made of sterile neutrinos with a mass in the keV 
range particularly interesting: 

(i) It is the lowest possible range of masses for fermionic 
DM. The Pauli exclusion principle applied to the DM parti- 
cles in the halos of dwarf spheroidal galaxies (such as Draco 
or Ursa Minor) imphes a lower bound on the particle mass 
Mom > 350 eV jTremaine & Gunn|[T979,) . Thus the sterile 
neutrino can be light enough to be a warm DM candidate (see 
below). 

(ii) Warm DM with a keV mass can alleviate the prob- 
lem of too many small subhalos inside the bigger dark mat- 
ter halos, and too sharp central density peaks in the galaxy- 
sized DM halos predicted in the Cold Dark Matter sce- 
nario (Bode et al."200II; lGoerdt et al.ll2006h . For example, the 
flat central radial mas s profi l e of the Fornax dwa rf spheroidal 
galaxy dGoerdt et all 120061; IStrigari et al] l2006h can be ex- 
plaine d if dark matter is w ar m with Mqm ^ 2 ke V. 

(iii^ lAsaka et alJ (l2005h and lAsaka & Shaposhnikov! (l2005h 
recently demonstrated that a simple extension of the Stan- 
dard Model by three singlet fermions with masses smaller 
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than the electroweak scale can accommodate the data on neu- 
trino masses and mixings, provides a candidate dark matter 
particle (in the form of the lightest sterile neutrino), and can 
explain the baryon asymmetry of the Universe. Such an ex- 
tension (dubbed i/MSM) can quantitatively explain these "be- 
yond the Standard Model" phenomena within a single consis- 
tent framework. It should be tested observationally, and one 
such test is the search for the sterile neutrino DM. 

(iv) For a DM particle with the mass in the keV range, one 
can obtain lower bounds on its mass by modeling the large 
scale structure formation. The power spectrum of the mat- 
ter density fluctuations derived from the Lyman-a for est data 
in the SDSS , spanning redshifts 2.2 < z < 4.2 (Sel jak et alJ 
l2006l;IViel et alj2006) constrains the mass of a warm DM par- 
ticle to the range > lOkeV (> 14keV in iSehak et alJ(l2006i) V 
However, a more conservative analysis using only the higher 
spectral resolution Lyman-a data and lower redshifts gives 
Mom >XkeW (Hansen et al. 2002; Viel et al. 2005). 

(v) Sterile neutrinos with a keV mass w ould have other in- 
teresting astrophysical applications (e.g.. iKusenkol l2006blal : 
iBiermann & Kusenkoll2006l;IStasielak et al.ll2006l) . 

T he sterile neutrino should possess a radiative d ecay chan- 
nel (iPal & WolfensteinI flQSl fearger et al.1 fTool (see €lT\ 
below). An emission line from neutrino decay has been 
searched for — so far unsuccessfully — in the X-ray 
spectra of various ty pes of astrophy sical objects , including 
cluster s of galaxies (lAbazajian et al.ll2001b: Bovarskv et al. 
2006d), the diffuse X-ray b ackground dDolgov & Hansen 



Way (|Boyarsky et al. 2006d e; Riemer-Sdren sen et al.l l2006). 
dwarf galaxies (Bovars kv et al.l I2006dllel) . and the M31 
galaxy (Watson et a l. 20061). Assuming that all of the dark 
matter is in the form of sterile neutrinos, nondetection of such 
af line places constraints on the mixing angle of the sterile 
neutrino as a function of mass in the range 1 keV < Mj < 
lOOkeV. To derive such constraints, one needs to know the 
mass of the DM in the field of view of the X-ray spectrom- 
eter, Mjj"^ (eq.|2] below). There are several ways of deriving 
these masses: 



2002tjBovarskv et al.' '2006a'), the DM hal o of the Milky 



Fig . 1. — Regions us ed for extracting the X-ray spectra overlaid on a Chandra X-ray 0.8^ keV image of IE 0657-56 (left panel) and its weak lensing mass 
map jClowe et alj|200d) (right panel). Tlie panel size is 12' X 12'. All point sources seen in the X-ray image are spatially excluded from the spectral analysis. 
The region that we refer to as SUB is the green circle centered on the western subcluster's mass peak, excluding the smaller black circle encompassing most of 
the X-ray gas bullet. The region PEAKS is a combination of SUB and the green circle centered on the eastern mass peak, again excluding the con'esponding X-ray 
peak (the bigger black circle). The region WHOLE is the big (r = 6') green circle, excluding only the gas bullet (the smaller black circle). 



(a) Modeling rotational curves of stars in galaxies or the 
velocity dispersion of galaxies in dynamically relaxed galaxy 
clusters; 

(b) reconstructing the density and temperature profiles of 
the hot intergalactic gas in relaxed galaxy clusters using X- 
ray data and determining the tot al mass from the assumption 
of hydrostatic equilibrium (e.g., 'Cavaliere & Fus c _ p-Femianol 
II976; Bahcall & Sarazin 1977; Vikhlinin et al. 200^^ 

(c) For sufficiently distant ( z > 0.1) clusters , gravitational 
lensing can be used CBartelmann & Schneider! 12001) . which 
does not require a cluster to be relaxed. 

All previous observational constraints were derived for 
nearby objects {z < 0.01), for which a combination of meth- 
ods (a) and (b) provided the DM density distribution andM^",^. 
It is important to cross-check these results using objects with 
masses determined by all methods, to minimize the system- 
atic uncertainties inherent in each method. 

In the present work, we use a distant object (z = 0.296, cor- 
responding to the luminosity distance Di=\ .530 Gpc for our 
adopted cosmology with h = 0.7, Vl^ = 0.3 and Q,/^ = 0.7), 
whose mass is determined via weak and stro ng gravitational 
lensing (IClowe et al.ll2006t iBradac et al.ll2006 ). This method 
gives the total projected mass within a given region in the sky, 
eliminating a number of uncertainties of other methods. It is 
the only mass measurement method applicable to this clus- 
ter, which undergoes a violent merger. The merger has also 
resulte d in a unique separation between the dark and visible 
matter (IClowe et al.| [2006), which we will utilize below. 

1.1. Radiative decay of sterile neutrinos 

A dark matter co mposed of sterile neutrinos should not be 
comple tely "dark" dPolgov & HansenI 120021; lAbazaiian et al.l 
1200 Ibl) . The sterile neutrino possesses a radiative decay chan- 
nel, decaying at a rate F into an active neutrino and a photon 
with energy E = Mjl (where Ms is the sterile neutrino mass). 
It is convenient to parameterize the interaction with the ac- 
tive neutrino in terms of the mixing angle sin^(26'). The decay 



rate F is then given by (iPal & Wolfensteinlll982HBarger et all 
[1991 : 

'^"^^;sin2(20)M5 



10247r4 ' 



= 1.38 X 10"^^sin^(26l) 
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For distant objects, the decay flux into a solid angle fifov 
(the spectrometer's field of view, FoV) is given by 



P _M'Z^ E 

^nivi — 



(2) 



where M|j°^ is the total mass of DM within this solid angle. If 
the object is at a redshift z, 
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Nondetection of any X-ray emission lines in the spectrum of 
a massive object that are not expected from its baryonic con- 
stituents (the hot intergalactic gas in the case of a galaxy clus- 
ter) can be used to place an upper limit on the flux from sterile 
neutrino decay. Eq. (O can then be used to constrain the pa- 
rameters Ms and 9. 

2. THE IE 0657-56 CLUSTER 

IE 0657-56 is an interesting object for constraining the 
brightness of the neutrino line, for several reasons. Its 
total mass is directly measured from gravitational lensing 
(IClowe et al.1 120061; iBradac et all l2006h. It also has a very 
long (450 ks) Chandra observation (lMarkevitchl2005h . which 
provides a high-statistic X-ray spectrum. While formal con- 
straints on the neutrino model that we will obtain below (® 
are not significantly better than those previously derived from 
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TABLE 1 

Masses within spectral extraction regions 



Region 


Total Mass," 


Gas mass, 
IO'^Mq 


DM mass, 
IO'^Mq 


area,* 
10-' sr 


Sub 


0.058 


0.007 


0.05 


1.00 


Peaks 


0.198 


0.034 


0.16 


3.55 


Whole 


1.46 


0.297 


1.16 


95.3 



" Masses from weak leasing 
* The solid angle of the region 



the ne arby X-ray clusters Coma and Virgo dBoyarsky et alj 
l2006d) . IE 0657-56 provides a significant improvement in re- 
liability, because of its directly measured total mass (whereas 
Coma and Virgo are both nearby unrelaxed systems, so their 
dark matter masses, and how much of it falls inside the instru- 
ment FoV, is uncertain). 

A unique feature of IE 0657-56 is a spatial separation be- 
tween the peaks of gas and dark matter density belonging 
to the two subclusters dClowe et alj l2006h . caused by their 
merger in the plane of the sky (Fig.[TJ. This enables us to try 
to exclude the spatial regions with the highest thermal X-ray 
contamination (and thus minimize the uncertainty of model- 
ing this component, see ^3. II below), while at the same time 
retaining the densest dark matter regions. Thus, we will use 
two regions in our X-ray analysis below, shown in Fig.[T] The 
region PEAKS takes advantage of the separation between the 
DM and gas and combines two circles centered on the mass 
peaks, excluding the two X-ray brightness peaks. We will 
also use a subregion of PEAKS that includes only the bullet 
subcluster mass peak, SUB, to illustrate the effects of uncer- 
tainties. The region WHOLE includes most of the cluster mass 
within r = 6' = 1.6 Mpc (a still bigger region will increase 
uncertainties of the mass and the detector background), ex- 
cluding only a small region at the X-ray brightness peak. 

To calculate the dark matter mass within the spectral ex- 
traction regions, we i ntegrate d the weak lensing map of the 
projected total mass (IClowe et al. 2006) and subtracted a rel- 
atively small contribution from the intracluster gas. We note 
that the mass near the cluster center derived from weak lens- 
ing is lower by a factor of about 2 compared to that derived 
by Bradac et al., who combined weak and strong lensing 
data (but whose fit is in fact dominated by the strong lens- 
ing data). Some of this discrepancy is expected, because the 
weak lensing approximation breaks down near the peaks of 
massive clusters that produce strong gravitational arcs (such 
as both subclusters of IE 0657-56); indeed, the peak densities 
in IClowe et al.l (i2006) are insufficient to produce arcs. Weak 
lensing is also insensitive to a dding a constant mass sheet 
dBartelmann & Schneiderl2001h . Strong lensing analysis may 
suffer from other types of uncertainties. The higher-mass 
iBradac etal . (2006) map is limited to the central r= 1.5' -2' 
region, insu fficient for our purposes, so we chose to use the 
IClowe et~al.i C2006.) mass map. This will result in conservative 
underestimates of the expected neutrino signal by a factor of 
up to 2. We will illustrate this uncertainty in the final results. 

The gas mass within our spectral regions is estimated from 
a three-dimensional model fit to the Chandra X-ray brightness 
and temperature maps. The X-ray emissivity of a hot gas (T 
8-20 ke V) in the Chandra energy band is determined mostly 
by the gas density, and depends weakly on temperature. Clus- 
ters are optically thin for X-rays, so in general, their gas den- 
sity can be reconstructed very reliably for symmetric clusters. 
The apparent axial symmetry of IE 065 7-56 allowed us to 
derive a gas model with a 10% accuracy iMarkevitch & et all 



( |2006|) . The gas contribution to the total mass in our spectral 
regions is about 25% (for the total mass from weak lensing), 
which we subtracted to obtain the dark matter masses given in 
Table[T] Given the small contribution of gas to the total mass, 
its uncertainty will be neglected. 

3. DATA ANALYSIS 

We use the deep Chandra (' Weisskopf et al."2002') ACIS-I 
observation of IE 06 57-56 described in Markevitch (2005); 
IMarkevitch & et al.l (|2006|) . The ACIS non-grating energy 
resolution is between 1 2%^% (half-power line width AE /E) 
in the £ = 1 - 8 keV range. The X-ray spectra, the blank- 
sky background spectra, and ACIS responses (that include 
all instrument effects such as mirror effective area, detector 
efficiency and energy resolution, to be applied to a model 
spectrum in order to compare it to the d ata) were derived as 
described in IMarkevitch & et all (l2006h and IVikhhnin et al.1 
(|2006|). Spectral fit ting was performed using the XSPEC pack- 
age (lArnaudlll996l) . The spectra for our three regions, SUB, 
PEAKS, and WHOLE, were well-fit (with reduced ~ 1 in 
all cases) with models cons isting of one or se veral thermal 
plasma components (APEC, ISmith et ani2001h representing 
the multi temperature intracluster gas. A useful property of 
thermal spectra is that a continuous range of gas temperatures 
produces a spectrum that can be adequately fit with a sum 
of just several discrete temperatures. Models for all regions 
were modified at low energies by the Galactic absorption with 
Nh = 4.6 X 10^" cm"^. An example of the fit for the region 
SUB is shown in Fig. |2l this fit has a reduced = 0.98 for 
152 d.o.f. 

None of the regions exhibits any emission lines other than 
those expected from hot gas (mostly the £ ~ 6.7 keV Fe line, 
redshifted to 5.2 keV). We use this fact to place constraints 
on the neutrino decay flux fo llowing the procedure p reviously 
used in this context by, e.g.. lBovarskv et al.l d2006a lHdllel). In 
particular, for each energy bin in the 0.8-9 keV range, we 
added a narrow emission line (a Gaussian line much narrower 
than the detector spectral resolution) to the thermal model, re- 
fit the spectra, and calculated a statistical upper limit on the 
line flux by increasing the line normalization until of the fit 
worsens by 9 (3(7 or 99.7% confidence level). To obtain con- 
servative upper limits, we allowed as much freedom for the 
parameters of thermal model as possible, including allowing 
the heavy element abundances (that produce the thermal emis- 
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Fig. 2. — Spectrum for the SUB region with the best-fit APEC model showin 
in red. For illustration, we show an additional narrow line at £" = 1 .3 keV (blue 
model line and residuals), which worsens the fit at a 3(T level. 
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Fig. 3. — Statistical upper limits (3cr) for the flux in a nonthermal, narrow 
emission line as a function of line energy, for our fitting regions. 

sion lines) to vary, thus letting the neutrino line mimic some of 
the thermal line flux at the respective energies. We note that 
the full number of counts in each bin (including the instru- 
mental background) is sufficiently high to ensure the Gaussian 
statistics, and is much higher than the resulti ng Umits on the 
line fl ux, so the use of Ax^ is appropriate (cf. lProtassov et all 
I2OO2I) . An example of an emission line that would correspond 
to a 3(7 upper limit is shown in Fig.|2] and the resulting statis- 
tical limits, in Fig. [3] 

3.1. Systematic uncertainties 

In addition to the statistical upper limits on the line flux, 
there are systematic uncertainties that has to be taken into 
account. First, the way how we n ormalize the ACIS back - 
ground using the high -energy band lMarkevitch et al.l (l2003h : 
iHickox & Mark evitchI (I2OO6) results in a 3% uncertainty of 
the normalization at the useful energies. Because the ACIS 
detector background has several prominent emission lines, 
such incorrect normalization may, for example, hide an emis- 
sion line coming from the sky, or create a spurions line. To 
take this into account, we varied the background normaliza- 
tion by 3% and repeated the fitting procedure. As expected, 
this leads to a noticeable increase of the allowed line flux only 
at rather high energies E >6 ke V, where the background in- 
tensity increases steeply. For the region WHOLE, limits with 
nominal background normalization and those of for normal- 
ization changed by ±3% are shown in Fig. H) These differ- 
ences were added in quadrature to the statistical limits at each 
energy. 

There is a more insidious uncertainty arising from the inac- 
curacies of calibration of the detector response and gain (the 
energy to spectral channel conversion). To assess this uncer- 
tainty, we have extracted a spectrum from the 880 ks ACIS 
observation of the the very bright Perseus cluster, and fit it us- 
ing the same calibration products (current as of summer 2006) 
and a model consisting of several thermal models as we use 
in this work, with all element abundances allowed to vary. 
The fit is shown in Fig. |5] Statistical errors in this dataset 
are mostly negligible. The fit shows systematic residuals at 
a 2-3% level of the model flux, some edge-like or even line- 
like, obviously caused by calibration inaccuracies (e.g., the 
feature around E = 2 keV is obviously due to a gain error). To 
take this uncertainty into account, we added 3% of the ther- 
mal model flux contained within the width of a Gaussian line, 
in quadrature to the statistical limits on the line flux. The re- 
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Fig. 4. — The effect of a ±3% systematic uncertainty in the ACIS back- 
ground normalization for the region WHOLE on our line limits. It is significant 
for large regions and at high energies; for smaller regions such as SUB and 
PEAKS, it is negligible (not shown). 

suit of adding these uncertainties for regions WHOLE and SUB 
is shown in Fig. |6] As expected, this uncertainty contributes 
mostly at lower energies, where thermal emission is bright 
(at high energies, the increasing statistical uncertainty starts 
do dominate). This is the uncertainty that can be minimized 
by observing "dark" matter clumps, such as our gas-stripped 
subcluster 

Finally, the biggest uncertainty, unrelated to the X-ray data, 
comes from the factor of 2 difference between the cluster 
masses determi ned from the weak and strong gravitat ional 
lensing analyses lClowe et al.l(l2006h : lBradac et alJ(l2006h . We 
will include it in the plot with results below (Fig.lTJ?). 

4. RESULTS 

Using eq. (O and masses from Table [T] (and, of course, 
the assumption that sterile neutrinos account for all of the 
dark matter), we convert upper limits on the neutrino line flux 
for our two regions into restrictions on sterile neutrino in the 
Ms — sm^(20) plane. They are shown in Fig.|7] The strongest 
constraint comes from the region WHOLE — the bigger mass 
of DM within the field of view turns out to offer a greater 
advantage than the reduced thermal contamination at the gas- 
stripped DM peaks. For the SUB and PEAKS regions, we plot 
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Fig. 5. — Chandra ACIS spectrum of the Perseus cluster from a 880 ks 
exposure, extracted from the 8' X 8' central region, excluding the very center 
(r < 1') with complex gas structure. The fit residuals illustrate the current 
calibration uncertainties. The residuals around 2-A keV are 2-3%. 
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Fig. 6. — The effect of a 3% calibration uncertainty in the model flux (that 
can take the .shape of spurious lines, see Fig. [3). It is significant for the line 
limits derived in regions with bright underlying X-ray emission (WHOLE), 
and negligible for regions avoiding the X-ray brightness peaks, such as (SUB). 

constraints for the conservative (perhaps excessively so) DM 
mass estimate from weak gravitational lensing (Fig.lTJj) along 
with a stronger estimate for the higher strong lensing mass 
(Fig. HI;). 

Although on average these results are about an order of 
magnitude weaker (in terms of sin^[2 01) than o ther recent 
limits (e.g., Boyarsky et al . 2006d; W at son et al.l 2006) ^ or 
more recent results ( Bo yarsky et al.l |2006e), they serve as 
an important cross-check. First, they are obtained f rom an 
object wit h z ~ 0.3, while previous res ults (Bovars ky et al.l 
l2006alHdl: lRiemer-S0rensen et al.ll2006t IWatson et al.ll2006l) 
were obtained for objects with z < 0.01 ^. Furthermore, the 
DM mass was determined via gravitational lensing, a method 
not applicable for nearby objects. This is important, as dif- 
ferent mass measurement methods are subject to different un- 
certainties, and using an object such as IE 0657-56 makes the 
constraints more robust. 

4.1. Dodelson-Wilson scenario 

Assuming that sterile neutrinos constitute all the DM 
and that they are produced in the early Universe via mix- 
ing with active neutrinos only, one should expect a rela- 
tion between the mass of sterile neutrino and its mixing 

' We note that the M3 1 constraints presented in IWatson et alj 120061) (as 
a straight line in the Ms -sin- (28) plane) should be very qualitative at M, > 
10-12 keV. They must worsen as our limits do, since the XMM effective area 
rapidly declines at the con'esponding energies, similarly to Chandra's. 

* When this work was in final preparation, a 



angle (iDodelson & Widrowl 11994 iDo lgov & HansenI 120021: 
lAbazajian et al.l !2001a: Abazajian 2006^! Combined with 
observational restrictions on the DM decay emission, this pro- 
vides an upper bound on the sterile neutrino mass. 

However, as sterile neutrinos do not thermalize in the early 
Universe, any such model relies on a number of assump- 
tions (including initial conditions at temperatures > 1 GeV 
and the absence of entropy dilution) dBovarskv et al.ll2b06dl; 
lAsaka et"ani2006ah . To define the boundary conditions, the 
knowledge of some "beyond the i/MSM " physics is needed. 
For e xample, it was shown in Shap oshnikov & TkachevI 
(l2006l) that all of the DM sterile neutrinos could have been 
produced by interaction with inflation. In such a scenario, 
the mixing angle can be arbitrarily small, even zero. This 
means that the sterile neutrino mass and mixing angle are 
not necessarily relate d, contra ry to the assu mption u sed in 
recent literature (e.g.. lAbazajia n 2006b; Wat son et al.l l2006). 
Upper limits on the neutrino mass placed in those works, 
based on the upper limits on the decay line flux, are not only 
model-dependent, but depend strongly on the initial condi- 
tions. Moreover, even assuming ad hock initial conditions 
implying that there was no DM at the temperatures > 1 GeV 
(which is hardly physically justified), the correct calcula- 
tion of the production rate requires calculation of the non- 
preturbative QCD contributions, which is still a subject of dis- 
cussion in the literature ( Asaka et al. 2006a b). 

Nevertheless, for the sake of comparison with other works, 
the intersection of our constraints in Fig. |7}j with the Mj - 
sin^(20) relation obtained in lAbazaiianl (l2006bh for the sim- 
plest DW model (with one sterile neutrino, assuming zero 
initial conditions and a particular form of QCD contribution) 
corresponds to an upper limit M, < 6.3 keV. 
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preprint IRiemer-Sorensen et alj I200fl) appeared, in which a Chandra 
grating spectrum of the cluster A1835 (z ~ 0.25) is analyzed. That analysis 
strongly underestimated the effect of the cluster angular extent on energy 
resolution of the grating spectrum, so we do not consider it here. 
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